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Abstract. We have carried out a search for 18-cm OH megamaser (OHM) emission with the 
Green Bank Telescope. The targeted galaxies comprise a sample of 121 ULIRGs at 0.09 < z < 
1.5, making this the first large, systematic search for OHMs at z > 0.25. Nine new detections 
Or of OHMs are reported, all at redshifts z < 0.25. For the remainder of the galaxies, observations 
constrain the upper limit on OH emission; this rules out OHMs of moderate brightness (Loh > 
10 3 L0) for 26% of the sample, and extremely bright OHM emission (Loh > 10 4 LO) for 73% 
of the sample. Losses from RFI result in the OHM detection fraction being significantly lower 
than expected for galaxies with Lir > 10 12 L . The new OHM detections are used to calculate 
an updated OH luminosity function, with $ oc Lq^ 66 ; this slope is in agreement with previous 
r \ , results. Non-detections of OHMs in the COSMOS field constrain the predicted sky density of 
OHMs; the results are consistent with a galaxy merger rate evolving as (1 + z) m . where m ^ 6. 
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c/2 ' 1. Introduction 

a • 

i— —i 1 OH megamasers (OHMs) trace of some of the most extreme physical conditions in the 
universe - in particular, the presence of an OHM signals specific stages in the merger 
process of gas-rich galaxies. OHMs can thus be used as probes of their environments, 
qq , both directly and indirectly. Characteristics of the maser emission itself can be used to 
measure extragalactic magnetic fields (via Zeeman splitting) and gas kinematics, while 
the presence of an OHM is a signpost for phenomena associated with galaxy mergers, 
including extreme star formation and merging black holes. OHMs are a unique tool in 
this respect due to their extreme luminosities and ability to be seen at cosmic distances. 

The total number of OHMs detected to date is still low. As of 2012, there are ^113 
OHMs published in the literature, with roughly 50% discovered in the Arecibo survey 
of Darling & Giovanelli (2000, 2001, 2002a). No OHMs at a distance of greater than 
1300 Mpc (z — 0.265) have been detected, and no large, systematic searches for high- z 
OHMs have been carried out. The association of OHMs with IR-bright merging galaxies, 
however, means that the density of OHMs is expected to be much higher at z ~ 1 — 2, 
coinciding with an increase in both merging rate and cosmic star formation. We have 
conducted a search for high-redshift OHMs using the Green Bank Telescope (GBT). 
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2. Sample selection and observations 

We constructed three samples of galaxies to search for high-redshift OHMs. None of 
the samples are fully complete or flux-limited, but draw on the catalogs of IR-luminous 
galaxies with well-defined redshifts that were available at the time. 
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Table 1. Properties of new OHM detections from the GBT survey 



Galaxy 


z 


apeak 
°1667 

[mJy] 


^^1667 

[MHz] 


\Ln\ 


IXnl 

L-^W J 


IRAS 00461-0728 


0.2427 


8.93 


0.937 


3.58 


2.94-3.30 


IRAS 01298-0744 


0.136181 


118.28 


0.785 


3.78 


3.10 


IRAS 01355-1814 


0.192 


6.99 


0.344 


2.81 


3.25 


IRAS 01569-2939 


0.1402 


40.67 


1.821 


3.93 


2.95 


IRAS 10597+5926 


0.196 


17.56 


0.915 


3.85 


3.13 


IRAS 12071-0444 


0.128355 


4.27 


0.598 


2.45 


3.07 


IRAS 16090-0139 


0.13358 


9.00 


3.110 


3.53 


3.35 


FF 0758+2851 


0.126 


4.67 


0.396 


2.67 


2.39 


FF 2216+0058 


0.212 


18.40 


0.347 


3.53 


2.58 



IRAS PSCz galaxies not visible from Arecibo: The first sample of galaxies consisted of 
IRAS sources included in the redshift catalog of the PSCz survey. Galaxies were selected 
according to similar criteria as in the flux-limited Arecibo sample, but included objects 
lying outside the declination limits of Arecibo (—1° < 5 < 38°). Galaxies selected for 
GBT observations had: a declination range of —40° < 6 < 0° or 6 > 37°, a redshift range 
of 0.10 < z < 0.25, and a lower-luminosity threshold of Le0fj,m > 10 1 L©. 153 galaxies 
in the PSCz met these criteria, of which 47 of the brightest candidates were observed 
according to the LST windows during the early commissioning phase of the GBT in 2002. 

Sub-mm and ULIRG galaxies from the field: The second sample of potential OHM hosts 
was assembled from flux-limited catalogs of ULIRGs at higher redshifts. We began with 
35 galaxies in the FSC-FIRST catalog (Stanford et al. 2000) which consists of targets 
detected in both the IRAS Faint Source Catalog and the 20-cm VLA FIRST survey. 
This was supplemented with 5 IR-bright galaxy pairs and 26 sub-millimetre galaxies. All 
galaxies have Lm > 10 n L©, with more than half having Lm > 10 12 L©. The highest 
redshift in this sample is at z = 1.55. 

Starburst galaxies from COSMOS: The third group of 19 OHM candidates was the 
last observed in our program, and made explicit use of the results from the first two 
samples. Targets were selected from the COSMOS field, a 2-deg 2 survey with deep spec- 



tral coverage from X-ray through radio wavelengths. Recent infrared (Willett et al. 2011" 



studies show that while OHMs are found in infrared-bright galaxies, the OHM fraction 
is much higher for starburst-dominated galaxies vs. AGN. Selection of OHM candidates 
began with COSMOS galaxies detected by Spitzer at 70 /zm, and then eliminating all 



targets except the LIRGs, ULIRGs, and HyLIRGs identified by Kartaltepe et al. (2010) 
We removed all galaxies identified as AGN or with Lm < 10 12 L©. Finally, we culled the 
target list based on the expected RFI conditions near the observed frequency bands. By 
limiting the observed OH frequencies to cleaner regions, we have a broader margin for 
error on the galaxy redshift. The two windows used are at v Q h s — 825 — 830 MHz and 
960 - 1005 MHz, equivalent to redshifted OH at 0.97 < z < 1.01 and 0.67 < z < 0.73. 

Observations: We observed the OHM candidates in several sessions at the GBT 
from 2002-2010, totaling approximately 150 hours. The majority of observations used 
the maximum available bandwidth of 50 MHz and 8192 channels. Integration times were 
selected with the goal of achieving ~ 1 mJy rms per channel for each galaxy. The data 
were reduced using standard routines in GBTIDL, including extensive flagging for RFI. 
After flagging, we fit the radio continuum around the expected line center with a poly- 
nomial function of order n — 5. This removed both intrinsic continuum structure from 
the target itself and any baseline structure not removed by the position-switching tech- 
nique. After stacking, the spectra were smoothed to a rest-frame velocity resolution of 
10 km s _1 . 



Searching for OHMs at z > 1 



121 



IRAS 00461-0728 




IRAS 01298-0744 

120 r 



IRAS 01355-1814 





Observing frequency [MHz] 

1600 1200 1000 800 700 



71600 72900 74200 
v npt [km/s] 



IRAS 16090-0139 



40300 40800 41300 
v., [km/s] 



FF 0758+2851 



57000 57500 58000 
v npl [km/s] 



FF 2216+0058 






38000 40000 42000 
v npl [km/s] 



IRAS 01569-2939 



37000 38000 39000 
v« [km/s] 



IRAS 10597+5926 



62500 63500 64500 
v., [km/s] 



IRAS 12071-0444 






41300 42100 42900 
v, [km/s] 



58000 59250 60500 
V, [km/s] 



37300 38150 39000 
v, [km/s] 







5.0 mjy : 


- 




— i i 

^^^^ 1 i 4. - : 


; J 










i. — — : 


= .A 




p J^— — * -Ta2 mjy ; 


- v "/ir "iJ 












if , , , 





3. Oo detection 
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Figure 1. Left: 18-cm spectra of the nine new OHMs discovered with the Green Bank Telescope. 
Right: Integrated OH luminosities all known OHMs, including data from the literature (red) 
and the new GBT detections (black). Upper limits from the GBT survey are indicated with the 
green arrows. The black curves show detection thresholds for a Av — 150 km s _1 OHM at a 
range of telescope rms sensitivities as a function of redshift. 

3. Survey results 

OH megamasers: Out of 128 galaxies observed for OH, we detected new OH megamaser 
emission in nine objects (Figure [IJ. Seven of the detections were PSCz galaxies from the 
first sample, while the other two were from the FSC-FIRST catalog in the second sample. 
All nine OHMs have redshifts near the lower end of the sample distribution, with the 
most distant lying at a redshift of z = 0.2427. The observed frequencies are in the range 
of 1300 — 1500 MHz, which is covered by the L-band receiver and has relatively little RFI 
compared to the GBT prime focus bands. The OH emission has been confirmed with 
GBT follow-up observations for five of the galaxies. We also confirmed the detection of 



the previously-discovered OHM IRAS 09539+0857 |Darling fc Giovanelli (2001) 

Table [T] lists the 18-cm radio properties of the new OHM detections. We give the 
galaxy's optical redshift (zhei), peak flux density of the OHM (Sfgg 7 ), the ratio of the 
integrated 1667 MHz emission to its peak flux density (Wiqqj), measured OH luminosity 
(log Loh), and predicted OH luminosity (L r ^f) based on the Loh — Lfir relationship 
in 



Darling fe Giovanelli (2002a) 



OH non-detections: 112 galaxies showed no confirmed detections of OH. The Lq" for 
each galaxy was (conservatively) derived assuming a boxcar line profile with a linewidth 
At; = 150 km s~ x and a 1.5cr detection. The rms was measured from baseline-subtracted 
continuum centered on the optical redshift of the galaxy and in a frequency range suffi- 
cient to cover the uncertainty in the optical redshift {Av b s = Az x v rest /[l + z]). 



4. Updating the OH luminosity function 

One of the goals of performing a search for OHMs at higher redshifts was to improve the 



measurements of the OH megamaser luminosity function (LF). Darling & Giovanelli (2002b) 



used the results of the flux-limited Arecibo survey to construct a well-sampled LF between 
10 2 2 Lq < Lqh < 10 3 ' 8 !/©, which followed a power law in integrated line luminosity of 
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4>[X] oc Lq° h 64 Mpc™ 3 dex^ 1 . This measurement was limited to a narrow redshift range, 
spanning 0.1 < z < 0.23. We constructed a new OH LF by combining the GBT and 
Arecibo OHM detections, using the 1/V a method and combining limits on both spectral 
line and continuum emission. We fit a power-law to all bins with more than one detection, 
yielding: 



log = (-0.66 ± 0.14) log Loh 
where $ is measured in Mpc~ 3 dex" 1 and Lqh in 



(4.91 ±0.41), (4.1) 

The new detections only change 
the slope measured by Darling & Giovanelli (2002b) by —0.02 and the offset by +0.10. 
Both values are well within the uncertainties of the combined LF, as well as that of the 
original Arecibo LF. 



Assuming the Malmquist-corrected relationship of Lqh oc L}^ from Darling fc Giovanelli (2002a) | 



this gives <&[Lm] oc L 



-0.83±0.18 
IR 



We compare this to the LF of ULIRGs in the local Uni- 



verse using the AKARI measurements of Goto et al. (2011) Folding in the OHM fraction 
derived from the combined Arecibo and GBT samples results in a slope of (—0.6 ± 0.2), 
a much shallower value than that measured from the AKARI galaxies (— 2.6 ± 0.1). This 
inconsistency may suggest either that OHMs are highly saturated or that the maser 
strength is only weakly correlated with global properties such as Ljji (Darling & Gio- 
vanelli 2002a). The OHM LF of |Briggs (2000) | assumed a quadratic OH-IR relation cor- 
responding to unsaturated masing; a decrease in saturation could potentially steepen the 
OHM LF up to a slope of -1.5. 



5. Constraining the evolution of the cosmic merger rate 

One of the ultimate goals of high-redshift OHM surveys is to use megamasers as tracers 
of the populations of merging galaxies as a function of redshift. Models of the merger 
rate as a function of redshift are typically parameterized with an evolutionary factor 
of (1 + z) m . The value of m, however, is not well-constrained {e.g., Kim & Sanders 
1998; Bridge et al. 2010). Sufficiently deep surveys of OHMs can provide an independent 
measurement of the parametrization of the merging rate. We calculated the predicted 
sky density of OHMs as a function of redshift: 



dN r , cD 



dVLdv v H v oy /(l + z) 3 n M + Q A \a lnlO 



.max )° — (LoH,min) a ) , (5.1) 



where a and b are parameters of the OHM LF from Equation 14.11 ($[Loh] — bL OH ), 
LoH,min is the minimum OH luminosity that can be observed at a given sensitivity 
level, and LoH,max is the upper physical limit on OHM luminosity (Darling & Giovanelli 
2002b). 

Using the upper limit of zero OHMs detected in the COSMOS field, we place an upper 
limit on the merger rate of m < 6. While this is still within uncertainties for the highest 
estimated values of m, the COSMOS limit is an important first step in using OHMs as 
an independent tracer. Further measurements of OH deep fields at higher redshifts will 
be crucial for a more accurate constraint. 



6. Future OHM searches 

The final results for the GBT OHM survey yielded a much lower detection rate (9/121 
= 7%) of megamasers than expected. Based on the high detection fraction for galaxies 
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with Lir > 10 12 Lq from the Arecibo survey, in addition to our careful selection of 
starburst-dominated galaxies, we had predicted a success rate of 20-30%. We attribute 
one of the primary causes of the low OHM fraction to be the sensitivity of the obser- 
vations. Figure [T] shows the upper limits for OHM candidates from our survey, along 
with the necessary rms sensitivity to detect OH lines as a function of rcdshift. To de- 
tect median-luminosity OHMs at z = 1 will require rms levels of 100 /xJy, and perhaps 
significantly more time devoted to individual targets. 

It must also be mentioned that RFI is a major culprit at v ^ s < 1 GHz, significantly 
restricting the redshift path and increasing the noise in each receiver band. Future OHM 
searches may benefit from interferometric observations (for which celestial RFI is uncor- 
rected from dish to dish), or from observations in more radio-quiet environments. The 
GMRT (India) and ASKAP (Australia) are promising instruments for interferometry, 
while the upcoming 64-m Sardinia Radio Telescope (Italy) and 500-m FAST (China) will 
be options to continue the search for OHMs at high redshifts. 
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